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Abstract 

We propose an explanation of the LSND evidence for electron antineutrino appearance based 
on neutrino decay. We introduce a heavy neutrino, which is produced in pion and muon decays 
because of a small mixing with muon neutrinos, and then decays into a scalar particle and a 
light neutrino, predominantly of the electron type. We require values of gm^ ~ few eV, g being 
the neutrino-scalar coupling and m^ the heavy neutrino mass, e.g. in the range from 1 keV 
to 1 MeV and g ~ 10 -6 -10 -3 . Performing a fit to the LSND data as well as all relevant null- 
result experiments, we show that all data can be explained within this decay scenario. In the 
minimal version of the decay model, we predict a signal in the upcoming MiniBooNE experiment 
corresponding to a transition probability of the same order as seen in LSND. In addition, we show 
that extending our model to two nearly degenerate heavy neutrinos it is possible to introduce CP 
violation in the decay, which can lead to a suppression of the signal in MiniBooNE running in 
the neutrino mode. We briefly discuss signals in future neutrino oscillation experiments, we show 
that our scenario is compatible with bounds from laboratory experiments, and we comment on 
implications in astrophysics and cosmology. 
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I. INTRODUCTION 



The LSND experiment at LANSCE in Los Alamos took data from 1993-1998 and ob- 
served an excess of v e PJ. This was interpreted as evidence of — > v e transitions. The 
KARMEN experiment at the spallation source ISIS in England was looking at the same 
appearance chanel in the years 1997-2001 at a slightly different baseline than LSND, but 
no positive signal was found j^f. Reconciling the evidence for Vpj— > ^appearance observed 
in LSND with the other evidence for neutrino oscillations p|, |4|, y, |y, |2J remains a challenge 
for neutrino phenomenology. It turns out that introducing a fourth sterile neutrino [8J1 does 
not lead to a satisfactory description of all data in terms of neutrino oscillations ' 19, 10] 
(see Ref. [ll| for a recent update) because of tight constraints from atmospheric U, so- 
lar 0, and null-result short-baseline (SBL) experiments j^, [12, 0, E3- I* 1 view of this, 
several alternative explanations have been proposed, some of them involving very specula- 
tive physics: four-neutrino oscillations plus neutrino decay three neutrinos and CPT 
violation a lepton number violating muon decay five-neutrino oscillations jl8| . 

four neutrinos and CPT violation |19j, a low-reheating-temperature Universe 0, CPT- 



violating quantum decoherenc e |2ll| . mass- varying neutrinos 22 1, and shortcuts of sterile 
neutrinos in extra dimensio ns [23j . A critical test of the LSND signal will come soon from 
the MiniBooNE experiment [2^, which is looking for — > v e appearance in a similar range 
of L/E u as LSND, L being the distance travelled by the neutrinos and E v their energy. 

In this work we assume the existence of heavy (mainly sterile) neutrinos, nh, with a 
small mixing with the muon neutrino. We denote their masses by mu- The LSND signal 
is explained through the decay of rih into a light neutrino, mixed predominantly with the 
electron neutrino. A natural way to introduce neutrino decays is by means of a term in 
the Lagrangian, which couples neutrinos to a light scalar, similar to the so-called Majoron 
models |2s| . This term can be related to the neutrino mass generation mechanism through 
the spontaneous breaking of a lepton number symmetry via a non-zero vacuum expectation 
value of the scalar field. Similarly to these models, we introduce an interaction term between 
rih, light neutrinos and a scalar. However, as we are only interested in the phenomenological 
consequences of such a term for the LSND experiment as well as for the other relevant neu- 
trino experiments, we assume that this term arises at low energy as an effective interaction, 
not necessarily related to neutrino mass generation. We show that in the simplest case of 
one heavy sterile neutrino, n 4 , for gm 4 ~ 1-10 eV, e.g. neutrino masses around 100 keV and 
neutrino-scalar couplings g ~ 10~ 5 , LSND data can be explained in complete agreement 
with the null-results of other SBL neutrino experiments. 

Light neutrino decay has been considered as an alternative solution to neutrino oscil- 
lations, to explain the atmospheric neutrino anomaly [3, [22} and the solar neutrino puz- 
zle 0, 01. Although the recent observation of an oscillatory behaviour in atmospheric 
neutrino j^J and KamLAND (3| data excludes a pure decay scenario as the possible ex- 
planation of these problems, a combined scenario of neutrino oscillations and decay is not 
excluded yet. We do not consider these possibilities here and we assume that atmospheric, 
solar and KamLAND neutrino transitions are explained by neutrino oscillations. 

The decay of an exotic neutral massive particle, called karmino, produced in pion decays 
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was considered as a possible explanation for an anomaly seen in the first data set of the 
KARMEN experiment j^], which however, was not subsequently confirmed j^. This 
hypothetical particle would move non-relativistically {(3 ~ 0.02) from the source to the 
detector. If it were a massive neutrino, strong bounds on the mixing and lifetime could 
be placed from anomalous pion and muon decays j^jj (for a recent review of the bounds, 
see Ref. Differently from this case, where no specific decay scenario was adopted, we 

assume here a Majoron-like type of decay. Furthermore, for typical values of masses, the 
heavy neutrino is relativistic at the energies in the LSND and KARMEN experiments, and 
we require it to decay before reaching the detector, with a lifetime much smaller than the 
one of the karmino. 

An explanation of the LSND result invoking the interplay between neutrino oscillations 
and neutrino decay has been considered previously in Ref. [l5|. A fourth neutrino is intro- 
duced with a mass of a few eV, which decays into a massless Goldstone boson (the Majoron) 
and into light neutrinos. The signal in LSND is provided by mixing, similar to (3+1) four- 
neutrino oscillations, whereas the decay serves to circumvent the constraints from the CDHS 
experiment. This approach requires very large couplings g ~ O(10 3 ) and a mixing of v e with 
the heavy mass state, and hence it appears to be in conflict with various bounds on neutrino- 
scalar couplings (see Section IVTT1 for a discussion). In contrast, in our model, the signal in 
the LSND experiment is provided entirely by the decay, and no mixing of v e with the heavy 
neutrino is required. This allows us to invoke neutrino masses in the 100 keV range, which 
in turn permits a small coupling of g ~ 10~ 5 , in agreement with existing bounds. 

The outline of the paper is as follows. In Section |H] we present the minimal decay 
framework and deduce the probabilities, which are needed to calculate the event rates. In 
Section ITTT1 we present our analysis of LSND and KARMEN data within this scenario, and 
in Section IIVI we show that the decay framework is not in contradiction with any other 
null-result experiment by performing a combined analysis; we compare the quality of the fit 
to the one for oscillations. In Section IS we discuss the prediction for MiniBooNE and other 
future neutrino experiments within the minimal decay framework, whereas in Section IVII 
we show that, allowing for, at least, two sterile neutrinos and complex couplings, it is 
possible to obtain a CP-violating interference between oscillation and decay amplitudes, 
which would suppress the neutrino signal in MiniBooNE, but at the same time provide 
the antineutrino signal in LSND. In Section IVII1 we argue that our model is consistent 
with existing bounds from laboratory experiments, and we comment on implications for 
cosmology and astrophysics. In Section fVIIII we present our final remarks. 



II. DECAY FRAMEWORK AND TRANSITION PROBABILITIES 

Let us consider the general case of N Majorana neutrinos. We take the three light 
neutrinos responsible for solar and atmospheric oscillations to be ^1,2,3 with masses 717,^2,3 ~ 
1 eV, while the heavy neutrinos have masses to^.., 3> ^1,2,3- Furthermore we assume for 
the scalar mass mi )2 ,3 ^ <C m 4i5j , such that the three light neutrinos are stable. Hence 
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the terms in the Lagrangian relevant to the decay (with scalar couplings 1 ) are given, in the 
mass basis, by 

£ = -y j gtdViLn hR (j) + h..c. , (2.1) 

l,h 

where here and in the following the ranges for the indexes are I = 1, 2, 3 and h — 4, 5, 

In general the coupling matrix Qhi will be complex. 

In the case of Majorana particles, neutrinos are identical to antineutrinos. Weak interac- 
tions couple to left-handed (chiral) neutrinos and right-handed (chiral) antineutrinos while 
the propagation states are those of definite helicity. In the relativistic case, where helicity 
and chirality are approximately the same, one can identify neutrinos and antineutrinos with 
helicity states up to terms of order m/E u . It follows from Eq. (j2.1J) that not only the decay 
n — > v + (f) can occur, but also n — > v + (p is possible [34]. Both processes are suppressed by 
terms of order rrif l /E nh . For the former process the suppression is due to angular momentum 
conservation, while for the latter the chirality- flipping nature of the interaction, Eq. (|2.1j) . 
is responsible for the suppression. 2 

In the approximation m; ~ 0, « and rrih <§C £ % , the differential decay rate of an 
rih of helicity r with energy E nh into a v\ of helicity s with energy E n is 

^nl-*vf{E nh ) 1 ... |2 



dE Vl 16icE* h 



\M rs \\ (2.2) 



where the indices r, s take the values '— ' for neutrinos and '+' for antineutrinos, and the 
matrix element is given by 3 

For relativistic neutrinos, the limiting values of the final neutrino energy are < E Vl < E nh , 
and hence, the partial decay rate for n r h — ► vf + cf) for both cases, helicity-flipping (r ^ s) 
and helicity-conserving (r = s), is the same and given by |2J] 

r H = r, . = l^l 2m ^ (2 a) 



The total decay rate of n h is then given by Y h = 2 T 



1 The case of pseudo-scalar couplings gives exactly the same results in the relativistic approximation we will 
use. It can also be shown that for processes in which all of the involved stares are on-shell, a derivative 
interaction term can be rewritten in pseudo-scalar form |33j. 

2 Throughout this work we will assume Majorana neutrinos. Here we just remark that in the case of 
Dirac neutrinos the decay is analogous. However, the decay products of the helicity-flipping channel are 
unobservable, since they are the right-(left-)handed component of the neutrino (antineutrino) fields, which 
do not participate in weak interactions. 

3 According to angular momentum conservation the matrix element Eq. H2.3(l can be written also as 
| A^r-s | 2 = \ghi\ 2 m\{l =F cos6*)/2 for r = s (r ^ s), where 6 is the angle between the spin of rih and 
the direction of vi in the rest frame of n^. 
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The flavour neutrinos v a [a = e, /i, r, si, S2, ■ ■ •) are related to the massive neutrinos in 
Eq. (HH) by 

"a = E UalUl + E ^o* 71 * • ( 2 - 5 ) 

The differential probability that an (anti)neutrino of flavour a with energy E Va is converted 
into an (anti)neutrino of flavour (3 with energy in the interval [E V0 , E Uf3 + dE U0 ] is: 



dP v r^ v s(E Va ) 



dE, 



mfL 
2E~ 



+ 



/ , U /3h u , 



( s )rr( r )* pyn 



+ W r .{E Ua ,E; 



"0/ 



<LL'- 



■ m l L 

'2E„_ 

dP% 



S(E Ua - E, 



up J &rs 



exp 



dec 



dL> 



b{E Va — E Vf} ) 5 rs 



(2.6) 

and = X* for r = +. In 



where we have introduced the notation = X for r = 
Eq. (|2.6|) , we have used the fact that the light neutrinos ui are not coherent with because of 
the large mass difference, while both v\ and are assumed to be coherent among themselves. 
The first term describes the oscillations of light neutrinos, whereas the second one describes 
oscillations of the heavy neutrinos weighted by the probability that they do not decay. These 
oscillation terms are only present in the helicity-conserving channel (r = s). The third term 
describes the appearance of decay products, and is present for the helicity-conserving as 
well as helicity-flipping channels. Here W rs (E Va , E Vj3 ) is the normalised energy distribution 
of the decay products: 



W r8 {E Va ,E v 



1 dT n r v b{E v 



hi 



dE, 



2®(E„ a 



E Vp ) x 



Eu p /E 2 Ua r = s, 

{E Va -E V0 )/El r^s, 



(2.7) 



where is the partial decay rate for —>■ v\ given in Eq. (J2.4j) . The relativistic approx- 
imation mi ~ allows us to factor out W rs (E Ua , E Vf3 ) in Eq. (j2.fi j) . The term dP„?L, v s/dL 
is the differential probability that the heavy component of the neutrino v T a decays at the 
distance [L, L + dL\ and the decay product interacts as a v%. Adopting the effective operator 
method of Ref. |35f it can be calculated as 



dP: 



dec 



dL 



E^E* 7 ^ -P 



(r)* i(r) 



'2E V 



exp 



VkL 



(2.8) 



Here Am is an effective amplitude for the decay of nn into ui, similar to the "appearance 
operator" of Ref. j^, and it is given by 



Am — ghiAh 



hi 



\9hi\ 



(2.9) 



Note that in the limit where the light neutrino masses can be neglected in the decay kine- 
matics, the real quantity A h is independent of the index I. Since for SBL experiments we 
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have AmfyL/E C 1 we can neglect oscillations of the decay products to derive Eq. ([2.8|h 
but we include the possibility of oscillations of the heavy states before they decay. 
In the simplest case in which N = 4, Eq. (|2.8|) becomes 



dP' 



dec 



* -^ a ,\ 2 W\ 2 A\e~^\ (2.10) 



dL 

where we have defined the couplings in the flavour basis by 



i 



III. LSND AND KARMEN 



Let us now calculate the v r — > v s e appearance probability relevant to the LSND and 
KARMEN experiments. Unlike other explanations of the LSND result based on sterile 
neutrinos, our model does not require a mixing of the electron neutrino with the heavy mass 
state. Therefore, for simplicity we assume U e ^ = 0. This implies that the first two lines in 
Eq. (|2.fij) disappear and using Eq. (j2.1()j) we find 

dP ^E^ ] = W rs (E^E u J \ \U^R e (1 - e~^ L ) , (3.1) 
where we have defined the branching ratios by 

R a = 1 -^, s 2 = £m 2 = £ l^l 2 - ( 3 - 2 ) 

^ I a=e,n,T,s 

We have used T 4 = 2g 2 A\, and for simplicity we neglect the dependence of g a on the helicity 
indices, i.e. we neglect complex phases of U a i and g& (compare to Eq. (|2.11jl ). The factor 
1/2 in Eq. (J3.1|) accounts for the fact that half of the initial neutrinos decay into neutrinos 
and half of them into antineutrinos. 

In LSND, neutrinos are produced by the decay at rest (DAR) of pions n + — > n + + and 
by the subsequent decay of the muons /x + — > e + +u e +u^. In the detector, the appearance of v e 
is searched for by using the reaction v e + P — * e + + n. In case of oscillations, only the from 
the muon decay contribute to the signal. In the decay scenario, in addition, the v ^ from the 
primary pion decay will give a contribution to the v e signal as well. Note that the v e from the 
muon decay will not contribute to the signal because of our assumption {7 e4 = 0. The pion 
decay leads to a mono-energetic beam with an energy of E^} = 29.8 MeV, whereas the 
muon decay gives z/ M with a continuous spectrum, tfi^Ep^), rising up to E? 13 * = 52.8 MeV. 
The total number of neutrinos is in both cases <fio — J dEp (fi^Ep ). The expected number 
of P e events with neutrino energy in the interval [E Pe , E Pe + dEp e ] is given by 



dE Pi 



/ dEp^Ep^) 



dEp e J~ dEp e 



(3.3) 
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FIG. 1: L/Ey spectrum for LSND for decay and oscillations compared with the data given in Fig. 24 of 
Ref. Q| . The hatched histogram shows the contribution of the line from the primary pion decay in the 
decay scenario. 

where a{E Ve ) is the detection cross section and C is an overall constant containing the 
number of target particles, efficiencies and geometrical factors. To calculate the observed 
number of events, we fold Eq. ()3.3|) with the energy resolution of the detector and integrate 
over the relevant energy interval. 

For the statistical analysis of LSND data we use the total number of signal events implied 
by the transition probability P = (0.264 ± 0.040)% 0. We have chosen an error such that 
we can reproduce Fig. 6 of Ref. |26j, which is based on an event-by-event analysis of the 
DAR data. Furthermore we include spectral data in the form of 11 data points given in 
Fig. 24 of Ref. (jj, where the beam excess is shown as a function of L/E v . We fit these 
data with a free normalisation in order to avoid a double counting of the information given 
already by the total number of events. In our analysis we do not use the data from decay in 
flight (DIF) neutrinos, since the appearance signal in these data is less significant than in 
the DAR sample. Moreover, since the total number of DIF neutrinos is at most two orders 
of magnitude less than the flux from DAR, one does not expect a relevant contribution 
to the v e appearance signal from DIF due to the helicity-flipping decay. 

Our analysis of LSND data gives a best fit value of Xmm = 5.6/9 dof for oscillations and 
Xmin = 10.8/9 dof for decay. The reason for the slightly worse fit for decay is the spectral 
distortion implied by the energy distribution of the decay products given in Eq. ([2.7)1 . In 
Fig. [U we compare the best fit spectra for oscillations and decay with the data. In the decay 
scenario we predict slightly too many events at low energies and too few events in the high 
energy part of the spectrum. Note that the contribution of the z/ M line from the primary pion 
decay is rather small. This follows from the fact that in the helicity-flipping decay mode, 
most of the decay products have a very small energy (see Eq. (|2.7jl ). which implies that they 
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oscillations decay 




sin 2 26 \U A \ 2 R 

FIG. 2: Allowed regions for LSND (solid) and KARMEN (dashed) at 99% CL, and LSND and KARMEN 
combined (shaded regions) at 90% and 99% CL. The left panel corresponds to neutrino oscillations and the 
right panel to the decay scenario presented in this work. 

fall below the detection threshold. 4 Although the fit in the decay scenario is slightly worse 
than for oscillations, the overall goodness of fit (GOF) is still acceptable. The \ 2 value cited 
above implies a GOF of 29%. Note, however, that a more refined analysis of the LSND 
spectral data might allow a significant discrimination between oscillations and decay. In the 
following we will assume that LSND data can be explained by our decay model and we shall 
proceed with the combination with other SBL data. 

First we discuss the compatibility of LSND and KARMEN in the decay scenario (see 
Ref. ji^ for a detailed analysis in the case of oscillations). Much as in LSND, also in 
KARMEN neutrinos are produced by the decay of pions at rest and the subsequent muon 
decay. However, in KARMEN detailed time information is available from the pulse structure 
of the proton beam and the excellent time resolution of the detector. Taking into account 
the muon lifetime of 2.2 /xs compared with the pion lifetime of 26 ns, a possible contribution 
of neutrinos from the pion decay was suppressed in the — > u e oscillation analysis by an 
appropriate time cut. Hence, in contrast to LSND the helicity-flipping decay of neutrinos 
from the line does not contribute to the v e signal in KARMEN, and the first term in 
Eq. (jH.Hjl is absent. For the statistical analysis of KARMEN data we are using the 9 data 
points as well as the expected background for the prompt energy given in Fig. lib of Ref. 0. 



4 Note that here we neglected the dependence of the couplings g a on the helicity indices due to complex 
phases in U a i and g±i according to Eq. If this additional freedom is taken into account, different 

branching ratios R e for the helicity-conserving and helicity-flipping decays can be obtained, which will 
modify the relative size of the line contribution. 
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In Fig. |21 we show the allowed regions for LSND, KARMEN, and the combination of both 
experiments for the case of oscillations and for the decay scenario. The comparison of the 
left and right panels of Fig. El shows that the compatibility of LSND and KARMEN in the 
decay scenario is at the same level as for oscillations. The reason is that both the oscillation 
phase, and the exponent r 4 L in Eq. Ijfl.lj) . have the same dependence on L/E u . This allows 
us to accommodate both the positive signal in the LSND experiment and the KARMEN 
null result, by taking into account the different baselines (L LSND = 35 m, Lkarmen = 18 m). 
In addition, as discussed above, the line gives a small contribution to the LSND signal, 
which is absent in KARMEN because of the time cut. 

A powerful tool to evaluate the compatibility of different data sets is the parameter 
goodness of fit (PG) criterion discussed in Ref. |3j|. It is based on the x 2 function 

XPG Xtot,min ,min ) 

(3.4) 

i 

where Xtot mm i s the x 2 minimum of all data sets combined and \ 2 min i s the minimum of the 
data set i. Applying this method to the case of LSND and KARMEN we find Xpg = 5-02 
for oscillations and Xpg = 4-97 for decay. These x 2 numbers have to be evaluated for 2 dof, 
corresponding to the two parameters in common to the two data sets (sin 2 29 and Am 2 for 
oscillations and \U^\ 2 R e and gm 4 for decay). This yields a PG of 8.1% for oscillations and 
8.3% for decay. 



IV. COMBINED ANALYSIS OF LSND AND NULL-RESULT EXPERIMENTS 

Let us now discuss the survival probabilities relevant to the analysis of SBL disappearance 
experiments. Under the assumption U e4: = in Eq. (|2.6|) . no v e disappearance is expected in 
SBL reactor experiments: Pp e ^p e = 1. However, since we need mixing of with the heavy 
states in order to explain the LSND signal, one expects some effect in disappearance 
experiments. To simplify the analysis, we will neglect in the following the appearance of 
in the decay products, i.e. we adopt the choice i? e « 1 and i? M)T ~ for the branching ratios 
of the decay. Then Eq. ()2.fi)l yields for the SBL survival probability 

Pf^t = C 1 - IM 2 ) 2 + IMV^ • (4-1) 

In the relevant L/E v range the most stringent bound on disappearance comes from 
the CDHS experiment [12], where the number of events in a near and far detector are 
compared. Therefore this experiment is only sensitive to the L dependent term in Eq. (JO}. 
Since this term enters only via |t/^4| 4 the sensitivity of CDHS to is rather poor in the 
decay scenario. Eq. (|4.1J) has to be compared to 

P?Z%? = 1 " 4|^4| 2 (l - IM 2 ) sm 2 (4.2) 

for oscillations in a (3+1) four-neutrino scheme. In this case L-dependent effects enter 
already at order jf/^] 2 , which leads to significantly stronger constraints than in the case of 
decay. Technical details of our CDHS data analysis can be found in Ref. 
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As shown in Ref. |39j , also atmospheric neutrino data provide a non-trivial bound on the 
mixing of with heavy mass states, i.e. on |f/ M 4| 2 . Following Ref. j^, in the case of (3+1) 
oscillations, the survival probability of atmospheric is given to a good approximation by 

^r^+iM 4 , (4-3) 

where Pi is an effective survival probability involving only oscillations of the light mass 
states mi, which is obtained by numerically solving the evolution equation within the Earth 
matter. Note that the parameter (tA^I 2 enters also in Pi. The recent update of a four- 
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gives 



neutrino analysis of atmospheric |3j , and K2K [5|] neutrino data performed in Ref. 
the bound 

\U^\ 2 < 0.065 (99% C.L.) , (4.4) 

which holds for m 4 ^> a/ Am^ TM ~ 0.05 eV, where Am ATM is the mass-squared difference re- 
sposible for atmospheric neutrino oscillations. In the decay scenario the survival probability 
of atmospheric neutrinos is modified to 

pATM.dec ^ pi + \U^ e -T,L ^ fl ^ (45) 

where in the last step we have used the fact that for decay rates relevant to LSND, we have 
T 4 L 3> 1 for atmospheric baseline and energy ranges. Comparing Eqs. (|4.3jl and ()4.5|) we 
find that neglecting the small term (t^l 4 the constraint obtained for oscillations applies 
also in the case of decay. Note that the decay will give a small contribution (suppressed 
by |t/^4| 2 ) to the P u r^ u s transition probability for atmospheric neutrinos, which potentially 
affects e-like events. In addition, in general one has to take into account also oscillations of 
the decay products with Am^, for atmospheric neutrinos, and some of the v T e produced in 
the decay will oscillate back to v r . However, under the assumption R e m 1, i? M)T « for the 
branching ratios, oscillations of the decay products with Am\ TM will be doubly-suppressed 
by \U^U e ^\ 2 . In the following decay analysis we will neglect such subleading effects in 
atmospheric neutrinos, and we will use as an approximation the x 2 (|^4| 2 ) fr° m Fig- 19 of 



Ref. [Ill] , which has been obtained for oscillations. A detailed investigation of atmospheric 
neutrino data within the decay scenario is beyond the scope of the present work. 

Now we turn to the global analysis of all SBL data, and we investigate the compatibility 
of LSND and null-result experiments within the decay scenario. In Fig. |3] we show the 
allowed regions from the appearance experiments LSND+KARMEN compared with the 
bound implied from disappearance experiments. In this bound we include data from the 
v e — > v e reactor experiments Bugey CHOOZ and Palo Verde 4l||, as well as 
data from CDHS 12], atmospheric neutrinos and the K2K long-baseline experiment (if, 
as described above. 

In the left panel of Fig. |3] we reproduce the well-known result (see e.g. Refs. jHl,|3^|) that, 
for (3+1) mass schemes, the effective oscillation amplitude sin 2 29 e s = 4|[/ e 4| 2 |[/ /J 4| 2 is tightly 
constrained from the disappearance experiments by the quadratic appearance of the small 
parameters \U e ^\ 2 and (f/^l 2 . In contrast, in case of decay, appearance and disappearance 
experiments are in perfect agreement, as is clear from the right panel of Fig. El there is a 
large overlap of the allowed regions for both data sets. The only relevant constraint from 
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FIG. 3: Allowed regions for LSND+KARMEN (solid) and SBL disappearance+atmospheric neutrino 
experiments (dashed) at 99% CL, and the combination of these data (shaded regions) at 90% and 99% CL. 
The left panel corresponds to neutrino oscillations in the (3+1) mass scheme and the right panel to the decay 
scenario presented in this work. The dash-dotted curve in the right panel shows the 99% CL constraint from 
CDHS. 

disappearance data comes from the bound from atmospheric neutrinos shown in Eq. (J4.4j) . 
As mentioned above no constraint arises from reactor experiments, and in agreement with 
the discussion related to Eq. (|4.1j) one can see in Fig. El that the constraint from CDHS is 
too weak to contribute significantly within the decay framework. The shaded regions in the 
figure show the allowed regions obtained from combining all data, where the total number 
data points is 

H(LSND) + 9(KARMEN) + 15(CDHS) + 60(B U gey) + l(CHOOZ) + l(Palo Verde) + 1 (ATM) = 98 • (4-6) 

The global analysis gives the following best fit parameters for the decay scenario 

\U^\ 2 = 0.016, gm 4 = 3AeV. (4.7) 

Hence, for neutrino masses in the range 100 keV coupling constants of order 10~ 5 are suf- 
ficient to make the decay fast enough. Let us quantify the quality of the fit further using 
three different statistical tests. 

1. First we compare the fit for oscillations and decay by using the absolute values of the 
global x 2 minimum. We find 

Xmin,osc 96.9 , Xmin,dec 88.3 . (4-8) 

Although the GOF implied by these numbers is good in both cases, thanks to the large 
number of data points (see Eq. (|4.6jl ). the A% 2 = 8.6 indicates that decay provides a 
significantly better description of the global data. 
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2. Next we use the PG [37] to test the compatibility of LSND with the rest of the data. 
Applying Eq. (|3.4|) to this case we find 



osc: Xpg = 21 - 8 > PG = 1.8 x 1(T 3 % 
dec: Xp G = 6.2, PG = 4.6 % 



(LSND vs rest) 



(4.9) 



where the x 2 values have been evaluated for 2 dof corresponding to the two parameters 
in common to the two data sets. The PG numbers show that for oscillations there is a 
severe disagreement between LSND and all the other experiments, whereas the PG is 
acceptable within the decay framework. The reason for the rather small PG even for 
decay comes from the slight conflict between LSND and KARMEN, which is present in 
any of the scenarios. Note that for five-neutrino oscillations in a (3+2) mass scheme, 
a similar analysis performed in Ref. |l8| yielded PG(3 + 2) = 2.1%, slightly worse than 
the one for decay. 

3. The better agreement of the data for decay becomes even more transparent if we 
test the compatibility of appearance and disappearance experiments, i.e. similar to 
Fig. we divide the global data into LSND+KARMEN and all the rest. Then the 
PG analysis gives 



In this analysis the conflict between LSND and KARMEN is removed, since they are 
added up to one single data set. Therefore, the above numbers confirm our conclusion 
from Fig. El that appearance and disappearance experiments are in perfect agreement 
in the decay scenario. 

Before concluding this Section, let us mention again that we have performed the analysis 
by setting {7 e4 = and R e = 1, = 0. This reduces the number of parameters, and 
hence simplifies the analysis considerably. Although we do not expect a significant change 
of our results by relaxing these assumptions, we stress that on general grounds this can only 
improve the fit of the decay scenario, since more parameters become available to describe 
the data. 

V. PREDICTIONS FOR MINIBOONE AND OTHER FUTURE EXPERIMENTS 
A. MiniBooNE 

A critical test of the LSND signal will come from the MiniBooNE experiment > which 
is currently taking data. This experiment looks for u e appearance in a beam of neutrinos 
with a mean energy of ~ 700 MeV at a baseline of L ~ 540 m. In the minimal decay 
framework discussed in the previous Sections we predict for MiniBooNE the appearance of 
v e events with a probability of the same order as in LSND, similar to the case of oscillations. 
Since the MiniBooNE detector cannot distinguish between neutrinos and antineutrinos, the 



osc: Xpg = 16 - 6 , PG = 2 - 5 x 10 -2 % 
dec: Xpg = 12 > PG = 55 % 



(app. vs disapp.) 



(4.10) 
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FIG. 4: Energy spectrum predicted for the v e appearance signal in MiniBooNE for decay with gm\ = 3.4 eV 
(shaded region) and for various values of Am 2 in the case of oscillations (solid curves). The blue/dark-shaded 
region shows the contribution of v e from the helicity changing decay. 



v e produced in the helicity-flipping decay n 4 — > v e (j) will also contribute to the signal. We 
calculate the differential number of events as 

oc / dE v J{E v ^) 

where a (a) is the total charged-current cross section for neutrinos (antineutrinos), and 
(j)(E u ) is the initial flux of z/ M , which we extract from Fig. 3 of Ref. [24]. 

The spectral shape predicted by Eq. ()5.1|) is shown in Fig. 0] Although the heavy neutrino 
mass state decays with equal probability into v e and v ei the antineutrinos give only a small 
contribution to the total signal. The reason is that, according to Eq. (j2.7J) . most of the decay 
products from the helicity-flipping decay will have low energies, for which the detection 
cross section is small. Moreover, the detection cross section for antineutrinos is roughly a 
factor two smaller than the one for neutrinos. Note that the spectral distortion implied by 
W rs (E Va1 E Vf3 ) will be less pronounced in MiniBooNE than in LSND because of the very 
different initial spectra. In the LSND experiment, the spectrum rises monotonically up to 
a maximum energy, whereas the MiniBooNE spectrum decreases with energy. Therefore, 
degrading the energy of the decay products has a smaller impact on the final spectrum in 
the case of MiniBooNE. 

In Fig. H] we show also the predicted spectral shape in the case of oscillations for several 
values of Am 2 . One observes that in principle the decay predicts a specific spectral shape of 
the signal, whereas the actual spectrum from oscillations depends strongly on Am 2 . Whether 
it will be possible to distinguish between the two models in the case of a positive signal in 
MiniBooNE depends on the available statistics, as well as on experimental factors such as 



dP„ 



XEu 



dE„ 



dP„ 



XEu, 



dE* 



o{E v 



(5.1) 
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the energy resolution for E v , which are not taken into account in Fig. 0] 

Another method to discriminate between decay and oscillations is provided by the disap- 
pearance channel in MiniBooNE. In the case of (3+1), as well as (3+2) oscillations sizeable 
disappearance is predicted, which should be observable in MiniBooNE j^]. In contrast, 
the signal is expected to be very small in the decay scenario. The disappearance search in 
MiniBooNE is limited to a shape analysis of the energy spectrum, since the normalisation 
of the total number of events suffers from large uncertainties 4^]. Therefore, the signal for 



decay is suppressed for the same reason as discussed above in relation with Eq. ()4.1|) and 
the CDHS experiment. 



B. Experiments looking for a non-zero value of U e 3 

In general the decay will contribute to the signal in experiments such as T2K or 



NOi^A 44J looking for a non-zero value of U e ^ by exploring the chanel — > u e . However, 
the appearance probability would be proportional to \U^\ 2 ~ 0.01, and thus at the edge of 
the sensitivity for those experiments. Therefore, it appears to be rather challenging first to 
observe an effect of our decay scenario, and second to disentangle it from oscillations induced 
by Uei- Much as in MiniBooNE, the signal from decay will have a characteristic spectral 
shape, as implied by the spectrum of the decay products. Moreover, in detectors capable 
of distinguishing neutrinos from antineutrinos the appearance of "wrong-helicity" neutrinos 
with low energies would be a generic prediction of the decay scenario. Nevertheless, this is 
also very challenging, because charge discrimination for electrons is a difficult experimental 
task. In this respect, an interesting possibility to explore the decay scenario might be to 
look for the appearance of v e from the beam of T2K in the KamLAND detector. Let 
us note that also the Minerz/a experiment Jiij], which will be placed ~ 1 km away from the 
target within the NuMI neutrino beam, or the K2K/T2K near detectors could be suitable 
places to look at these effects, although these beams might suffer from a large intrinsic v e 
(z/ e ) background. 

Another manifestation of the decay model could be the observation of — > v e appear- 
ance in accelerator experiments such as T2K or NOz/A, but no corresponding signal for v e 
disappearance in reactor experiments such as Double-Chooz This would imply very 

small values of U e z, such that v e disappearance is suppressed, whereas the — > u e ap- 
pearance signal is dominated by the decay. Note, however, that a positive signal for v e 
disappearance in reactor experiments cannot exclude the decay scenario, since our model is 
perfectly compatible with a sizeable value of 

In summary, we stress that a direct measurement of the decay effects is quite challenging 
for present and near-future experiments looking for a non-zero value of U e s, because of the 
rather small mixing, (t^] 2 , which is needed in order to accommodate the LSND result. 
Therefore, if the LSND signal should be confirmed by MiniBooNE a new experiment at a 
stopped pion neutrino source as proposed in Ref. [4jj could be an optimal experiment to 
confirm or exclude our decay model. 
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VI. CP- VIOLATING DECAYS AND THE SIGNAL IN MINIBOONE 



In this Section we extend our model and consider the case of two heavy neutrinos, i.e. N = 
5 massive neutrinos. We will show that in this case it is possible to obtain an interference 
between oscillation and decay amplitudes [35], which may lead to CP violation in the decay. 
In this way the signal in MiniBooNE (running in the neutrino mode) can be significantly 
suppressed with respect to the LSND signal from antineutrinos. 

Using Eq. (J2.8)) for N = 5 and performing the integral over L, one finds 

»l dP dcc 
/" l " r „ = / dl. 1 - 

U ah \ 2 R h p(l-e- T » L ) 



J dl. 
1 

2 



+ R/3 \U a 4\\U a5 \ cos7 cos(5 + 7)-e rL cos(5 + 7 + r osc L) , (6.1) 
with the following definitions for the branching ratios 

d _ n \9hp\ 2 A 2 h - _ - _ r 4 + r 5 

= * — p , Up = 2 = , 1 = , (b.2J 

1 h 1 ^ 

and phases 

5 = axg{U* A U a5 g^glp) , tan 7 = ^ , T osc = — ^1 . (6.3) 

The first term in Eq. (|6.1j) describes the incoherent decay of the two heavy mass states in 
analogy to Eq. (jSHJ), whereas the second term results from the interference of the amplitudes 
for decay and oscillations of the two heavy neutrino mass states and n*,. The phase 5 
is the analogue of the Dirac phase leading to CP violation in the standard three neutrino 
oscillation framework. In fact, in the case of antineutrinos, U a h and gh,p have to be replaced 
by U* h and g^g, which implies 5 — > —5, and we obtain the CP asymmetry: 



dec D<iec 

TL 



AP EE P- - - P 



2Rn \ U a A \U a $\ cos 7 sin 5 



sin7-e- iyj sin(7 + r osc L) . (6.4) 



One observes that necessary conditions to obtain CP violation are 5 7^ 0,7r and tan 7 ~ 1, 
i.e. T ~ r osc or ^m|/87r ~ Amj 4 . Using gh m h ~ few eV, the last condition shows that the 
heavy neutrinos (m^ > keV) have to be highly degenerate. 

In Fig. we show a numerical example for the Pp e< ^ Pe and P„ e< ^ Ue transition probabilities 
for an L/E u value relevant to LSND and MiniBooNE. Clearly, in this example the neutrino 
signal is strongly suppressed with respect to the antineutrino signal for Am^ ~ 1 eV 2 . 
This result confirms the previous observation that in order to obtain CP violation, the two 
heavy neutrinos have to be very degenerate, with masses in the keV range and a mass 
difference in the eV range. Whether this mechanism allows indeed to accommodate the 
LSND signal with a possible null-result of the MiniBooNE experiment in the neutrino mode 
has to be investigated after the release of the MiniBooNE data by a combined analysis of 
the two experiments within this framework. In that case, MiniBooNE data taking with an 
antineutrino beam would be necessary in order to test the CP-violating decay scenario. 
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FIG. 5: The transition probabilities Pp°£ +Pe and P^ e 1+ Ve according to Eq. Itj.l|) . The shaded region 
corresponds to the la range for P p _> Ce observed in LSND. The chosen parameter values are L/E u ~ 
0.75 m/MeV, |L/ m4 | 2 = W^ 2 = 0.003, R 4e = R 5 e = R e = I, T 4 = T 5 = T with g h m h = 3.4 eV, and 
6 = it/2. 

VII. CONSTRAINTS FROM LABORATORY, COSMOLOGY, AND ASTRO- 
PHYSICS 

Mixing between active and sterile neutrinos 5 , and couplings between active neutrinos and 
a light scalar, have been extensively studied, both in laboratory experiments and, for their 
implications in the evolution of the early Universe and of supernovae. For the discussion 
of these constraints we will focus on the minimal scenario necessary to explain the LSND 
signal: we require only mixing of with the heavy mass state, l?/^] 2 ~ 0.01, whereas we 
set £7 e 4 = U T 4 = 0. Furthermore, the assumption R e — 1, R PiT = for the branching ratios of 
the decay implies that only two elements of the (symmetric) coupling matrix in the flavour 
basis, g a/3 , are non-zero: g es ~ g and g efM ~ U^g. 

First we note that in our model solar neutrino oscillations are completely unaffected. 
Because of the assumption U ei = 0, no n h component can be produced in the Sun, and since 
the vi in our scenario are stable, there is no decay of solar neutrinos. For the very same 
reason the decay has also no effect in the KamLAND experiment, and these data are entirely 
explained by oscillations of the light neutrinos. On the other hand, effects in atmospheric 
neutrino experiments have been discussed in Section ITVl 



5 See Ref. |48| for a recent analysis of various bounds on sterile neutrino mixing. 
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A. Laboratory bounds 



In general, the mixing of a heavy neutrino leads to contributions to the effective neutrino 
masses in neutrinoless double-beta decay ji^ and tritium beta decay j5(J experiments. Note, 
however, that in our scenario there will be no effect of the heavy neutrino in such experiments 
because of our assumption U ei = 0. The coupling between v e and a light scalar would 
induce double-beta decay with the emission of one or two scalars, with a spectrum for the 
two electrons which is distinguishable from the one of the two-neutrino double-beta decay. 
Relatively strong bounds for single scalar emission apply only to the coupling g ee , which 
can be arbitrary small in our scenario. The couplings g ea for a/e contribute in principle 
to double-beta decay with the emission of two scalars. However, in this case the limits are 
very weak [^: 9 < 

The decay of pions and kaons has been used to set bounds on the mixing of heavy 
neutrinos (see Ref. ji^l and references therein). If a massive neutrino with m 1)2j 3 were 

produced in such decays, the energy spectrum of the muon would present an additional 
monochromatic line. No positive signal was found, leading to strong constraints on the 
mixing for neutrino with masses ^ 1 MeV. Mixing of neutrinos with smaller masses are 
compatible with these bounds. Furthermore, light scalar emission was not observed in pion 
and kaon decays 00. The most stringent bounds are of order g 2 < few x 10 5 , much 
too weak to constrain our model. 

B. Supernova bounds 

To estimate the effect of our model in thermal environments such as in a supernova and 
the early Universe let us compare the rate of processes induced by the Lagrangian Eq. (j2.1j) 
to the one of weak processes of the type ve — > ve. From dimensional considerations, it 
follows that cr wea k ~ GpE 2 , where E is the typical energy of the involved particles. Using 
cr(2 <->• 1) ~ g 2 m 2 h / E A and cr(2 <-> 2) ~ g A /E 2 , where 2^1 denotes processes like Vi<ft <-» n^, 
v i n h <P, and 2 <-» 2 indicates processes like V\V\ <-> <p<f), v\Vi <-> n^n^, we find 6 



In addition to these reactions induced by the scalar coupling, heavy mostly-sterile neutrinos 
are produced via oscillations due to their mixing with the active ones. 

If heavy neutrinos and light scalars were produced in the core of a supernova and could 
escape freely, they would carry away a large amount of energy, substantially modifying the 
supernova evolution. Together with the observations from supernova SN1987A, this energy- 
loss argument has been used to constrain sterile neutrino mixing J^J as well as Majoron 



6 For a careful analysis of these processes for active neutrinos see Ref. 




(7.1) 



(7.2) 
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coupling, see e.g. Refs. |33j, |55|. In our scenario the mass of the heavy neutrino is large 
compared with the matter potential within the supernova core. Therefore, we can treat 
mixing, as well as the reactions involving the scalar and rih, as in vacuum. From Eq. 1)7.1)) it 
follows that for typical energies in the supernova cores, E ~ 10 MeV, 2 <->• 1 reactions are 10 4 
times faster than weak interactions. This implies that rih and <fi are strongly coupled to the 
active neutrinos, and hence, they are trapped within the neutrinosphere, avoiding any energy 
loss due to particles escaping from the core. Let us add that there might be additional effects, 
such as modification of lepton number or neutrino flux distortions 0,0. The analysis of 
such effects requires detailed studies, which are beyond the scope of the present work. 



C. Cosmological bounds 

Constraints on the number of relativistic degrees of freedom present at the time of Big 
Bang Nucleosynthesis (BBN), usually parametrized by the number of neutrino flavours N u , 
are relevant when light particles such as sterile neutrinos are introduced |5^ |. Analyses 
of cosmological data, including light-element abundances lead to the limits 2.3 < N u < 
3.0 (95% CL) 0]. A recent analysis, that uses a new assessment of the primordial 4 He 
abundance and its uncertainty, has substantially relaxed the bounds on N u > 3: SN U < 1.44 
(95% CL) where 5N V = N v — 3. The discrepancy between these two exemplary results 
illustrates that such bounds have to be considered with care, because of the systematical 
uncertainties in light-element abundances, and the dependence on input priors. 

In the early Universe, rih of our decay model are initially generated by mixing with 
active neutrinos [fiiJ]. As soon as rih are produced, they are thermalized together with the 
scalar due to the very fast 2 <-> 1 reactions. Hence, one additional neutrino and one scalar 
degree of freedom will be present during BBN, leading to N v = 4.57, which appears to 
be disfavoured by the limits quoted above. 7 However, to draw reliable conclusions from 
BBN considerations a more detailed analysis is required, since there are many effects which 
might play a role. For example, in the presence of a small chemical potential for v e the 
limit on extra relativistic degrees of freedom becomes much weaker, N v < 6.5 jHgjl, and 
hence N v = 4.57 is allowed. Moreover, a detailed simulation of BBN in the presence of 
r-neutrino decay performed in Ref. showed that the effective N u depends in a rather 
non-trivial way on model parameters, and values of 5N U > as well as < may be obtained. 
Let us mention also the possibility to avoid the thermalization of rih and at BBN, by 
suppressing the initial production of rih through mixing. This could be achieved e.g. through 
a large lepton asymmetry of order 10~ 3 [62], a matter potential induced by neutrino-Majoron 
interactions 63], or by assuming a low reheating temperature [20]. The 2 <-> 2 processes 
would not be efficient in producing and thermalizing sterile neutrinos and scalars as far as 
they are slower than the expansion rate at BBN. Couplings g ^ few x 10~ 6 guarantee that 
these scatterings are out of equilibrium for T ^ 100 keV (see Eq. (|7.2jl ). 



7 Note that also in four- or five-neutrino oscillation schemes the sterile neutrinos are thermalized within the 
standard BBN scenario 1571. 
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As can be seen from Eqs. ()7.1|) and (|7.2|) the scalar interactions have a recoupling form, 
instead of the usual decoupling form of weak interactions, i.e. their strength grows as the 
Universe evolves, which may have some effects in later epochs of the Universe. For instance, 
if a light scalar couples to neutrinos before the recombination era, effects like a delayed 
matter-radiation equality and enhanced damping of acoustic oscillations at higher I can 
occur because the radiation energy density differs from that of the standard scenario. More- 
over, if one or more neutrinos are scattering during the eV era rather than free-streaming, 
a uniform shift of the CMB peaks to larger I would be manifest [64]. 

Finally, we note that our scenario, with the scalar heavier than the light neutrinos, does 
not give rise to a neutrinoless Universe which is disfavoured by structure formation 
arguments 



VIII. SUMMARY AND CONCLUSIONS 

We have presented an explanation for the LSND evidence for u e appearance based on 
the decay of a heavy sterile neutrino into a light scalar particle and light neutrinos. We 
assume a small mixing of the heavy neutrino mass eigenstate with of order jf/^] 2 ~ 
0.01, such that a small n.4 component is contained in the initial beam produced in the 
LSND experiment. On the way to the LSND detector the decays into the scalar and a 
superposition of light neutrino mass eigenstates. If this final state contains a large fraction 
of the v e flavour, it will give rise to the u e appearance signal observed in LSND. Taking into 
account the energy distribution of the decay products, one obtains a characteristic prediction 
for the spectral shape of the LSND signal, with more events at the lower part of the spectrum. 
Comparing this prediction with the spectral data of LSND, we find that the fit of the decay 
is slightly worse than for oscillations, although still acceptable (xlecay = 10.8/9 dof). A 
more detailed investigation of the LSND spectral data might allow to distinguish between 
decay and oscillations. 

We have performed a global fit to all relevant data, including null-result short-baseline 
reactor and accelerator oscillation experiments. The agreement of the — > v e appearance 
experiments LSND and KARMEN in the decay scenario is at the same level as in the 
case of oscillations. Much as the oscillation phase, also the decay exponential depends on 
L/E u , which allows us to reconcile the two results due to the different baselines. However, 
unlike the case of (3+1) four-neutrino schemes, our decay model is in complete agreement 
with the constraints from disappearance results. Using the so-called parameter goodness 
of fit (PG) to test the compatibility of appearance (LSND, KARMEN) and disappearance 
(Bugey, CHOOZ, CDHS, atmospheric neutrino) experiments, we find PG = 55% for decay 
but only PG = 0.025% for (3+1) oscillations. Testing the compatibility of LSND and all 
the null-result experiments, we find PG = 4.6% for decay, which is slightly better than the 
PG = 2.1% obtained in Ref. 3] for a (3+2) five-neutrino oscillation scenario. In addition, we 



have shown that our model is consistent with present laboratory bounds and have discussed 
implications for supernova evolution, BBN and CMB observations. 

In the minimal version of the decay model we predict a signal in the upcoming MiniBooNE 
experiment corresponding to a transition probability of the same order as seen in LSND, and 
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we have discussed characteristic signatures of our scenario, which may allow to distinguish 
it from oscillations in MiniBooNE or other future neutrino experiments. Furthermore, we 
have shown that within an extension of the minimal decay model it is possible to introduce 
CP violation in the decay through an interference term between decay and oscillations 
of two nearly degenerate heavy neutrinos, with masses in the 1 keV-1 MeV range and 
Am| 4 ~ 1 eV 2 . This can lead to a suppression of the signal in MiniBooNE running in the 
neutrino mode, and simultaneously accounting for the antineutrino signal in LSND. Hence, 
this model can only be excluded if MiniBooNE does not find a positive signal either in the 
neutrino mode or running with antineutrinos. 
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